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ABSTRACT 

We present a numerical study of the response of a thick accretion disc to a localized, 
external perturbation with the aim of exciting internal modes of oscillation. We find 
that the perturbations efficiently excite global modes recently identified as acoustic 
p-modes, and closely related to the epicyclic oscillations of test particles. The two 
strongest modes occur at eigenfrequencies which are in a 3:2 ratio. We have assumed 
a constant specific angular momentum distribution within the disc. Our models are in 
principle scale-free and can be used to simulate accretion tori around stellar or super 
massive black holes. 
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1 INTRODUCTION 

Accretion discs around compact objects (black holes and 
neutron stars) are natural laboratories for the study of 
strong gravitational fields and their effects (such as Lense- 
Thirring precession). From the physics that lies behind 
these accretion processes and the accompanying radia- 
tion, it is in principle possible to determine important 
parameters of the central object, such as mass and spin 
for the black hole cas e ([Abramowic^fc Khi zniakl l200lt 
IWaeoner. Silbergleit fc Ortega- RodrHmezT 1200 ll) . and to 
place constraints on the eq uation of state of dense matter in 
the c ase of neutron stars llKluzniak. Michelson fc Waeonerl 
Il990l) . 

It has also been pointed out that these objects 
can be potential sources of gravitational waves, aris- 
ing from differe nt processes, such a s nonaxisymmet- 
ric instabilities jZurek fc Gorskil Il989l) . nonaxisymmet- 
ric t ori powered by the spin energy of the black 
hole llva n Putter] |200l|) , magnetized thick accretion discs 
llMineshige. Hosokawa. Mashida fc Matsumotol [2002), and 
global oscil lations which induce a tim e-varying mass 
quadrupole JZanotti. Rezzolla fc Fonll2003l) . The energy re- 
leased due to accretion from such systems, if dynamically 
stable, has also been proposed as a mechanism for the pro- 
duction of cosmological gamma-ray bursts. In this case, 
super-Eddington accretion rates in flows cooled by neutrino 
emission feed the central black hole while powering a rela- 
tivistic outflow, which produces the GRB. Such a disc pos- 
sibly forms follo wing core collapse in a massive rotating star 
JWooslevlll993l) or after the merger of two compact objects 
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in a tight binary llNaravan. Paczvhski fc Piranl ll992). One 
mechanism that could potentially produce the collapse of 
such a disc onto the black ho le is the runaway radial insta- 
bility. This was identified bv lAbramowicz. Calvani fc NobiiH 
( 1983) and appears in discs with constant angular momen- 
tum distributions which overflow their Roche lobes at high 
rates, and thus increase the mass of the black hole sub- 
stantially. This alteration of the potential is a runaway pro- 
cess, and can destroy the entire disk within a few dynamical 
times. However, studies of tidal disr uption and torus forma- 
tion from compact binary mergers jRuffert fc Jankal 119991 : 
lLeell200ll) show that the distribution of angular momentum 
in the disc is far from b eing constant, and this can i n fact 
suppress the instability llDaigne fc Mochkovitcblll997l) . 

In the context of low-mas s X-ray binaries, obser- 
vations performed with RXTE llvan der Klisl koOO) have 
shown that there are millisecond oscillations in systems 
containing neutron stars or black holes surrounded by 
an accretion disc. An important result is that in at 
least four black hole sources (H1743-322, GRO J1655-40, 
XTE J1550-564, GRS915+105), two apparently stable 
peaks in the power spectrum appear at frequencies 
in the hH z range in a 3:2 ratio. A similar res ul t was 
found by lAbramowicz. Bulik. Bursa fc Kluzniakl j2003l) 
for the neutron star source Sco X-l. This lends sup- 
port to the resonance mode l originally proposed by 
lAbramowicz fc Kluzniakl (1200 J) , and further developed 
in terms of parametric resonance in a th i n dis c by 
Ajararnowicz jKaras. Kluzniak. Lee fc Rebuscol |2003) and 



R^btr^coll2004T) ~In this interpretation, the frequencies reflect 



epicyclic motion of perturbed flow lines in the accretion disc, 
or combin ations between these and a fixed, perturbation 
frequency llKluzniak. Abramowicz. Kato. Lee fc Stergioulasl 
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120041: iLee. Abramowicz fc Kluzniakl |2004), perhaps 
due to the stellar spin in neutron star sources 
jWii nands. van der Klis. Homan. Chakrabartv. Markwardt fc Morganl 
2003). Pressure coupling allows resonances to occur and 
manifest the mselves in the X-ray lightc u rve. 

Recently, IZanotti. Rezzolla fc Fontl i2003T) have shown 
that an extended torus can respond to external per- 
turbations in a global fashion, in a series of modes > 
whose frequencies follow the sequence 2:3:4:... . These 
are attributed to acoustic p-modes within the torus, ex- 
cited by an impulsive perturbation. Follow-up analyti- 
cal and numerical work ( Rezzolla, Yoshida fc Zanottill2003t 
iMontero. Rezzolla fc Yoshidall2004h has extended these re- 
sults, by calculating the set of corresponding eigenfrequen- 
cies (in height -integrated discs) and also by investigating the 
effects of different background metrics (e.g., Schwarzschild 
vs. Kerr for a rotating hole). In all cases where a numerical 
experiment was carried out, the perturbation was impulsive 
and global, affecting the entire torus. This idea has now been 
advanced as an explanation for the kHz QPOs in low-mass 
X-ray binaries containing the black hole candidates where 
the 3:2 frequency ratios have been reported, as ment ioned 
above iRezzolla. Yoshida. Maccarone fc ZanottSl l2003) 

In this Letter we show that it is possible to excite these 
modes in a thick accretion disc for a perturbation that is 
local, only affecting a small portion of the disc strongly. A 
frequency ratio of 3:2 for the two strongest modes is ap- 
parent. The lower frequency itself is related to the radial 
epicyclic value for test particles in circular motion, shifted *e 
to a lower frequency because of the finite extent of the torus. 
The second frequency maintains a 3:2 relation with the first. 



2 INITIAL CONDITIONS AND NUMERICAL 
METHOD 

2.1 Hydrostatic equilibrium for a thick torus 

We construct tori which are in hydrostatic equilibrium, have 
low mass, m < Mbh, and a radial extension L ~ R, where R 
is the distance separating the torus from the central mass, 
Mbh- They are thick in the sense that their vertical ex- 
tent, H, is comparable to L. We neglect the self-gravity of 
the torus, and additionally assume azimuthal symmetry. A 
polytropic equation of state, P = Kp 1 , has been used for the 
construction of initial conditions. Integrating the equations 
of hydrodynamics, it is possible to write: 

P 



= $ e + $0 = const, 



(1) 



7 

7-I/3 

where <I> e is the effective potential. Here $0 can be in- 
terpreted as a filling factor of the effective potential well. 
Through its variation tori of different sizes are constructed 
(Fig. 1 shows their cross sections over one half the r — z 
plane). The gravitational potential of the central mass 
is computed with the pseudo-Newtonian expression of 
IPaczvhski fc Wiital Jl980ft : 



—GMbh 
R-r a 



(2) 



which describes the behaviour of a test particle in a strong 
gravitational field and, importantly for our purposes, repro- 
duces the existence and positions of the marginally stable 




Figure 1. Surfaces of zero density and pressure, in meridional 
cross section, for various tori constructed according to the pro- 
cedure described in the text. Each curve corresponds to different 
values of the filling factor <E>o in units of c 2 /2. 



and marginally bound orbits in General Relativity (r g = 
2GMbh / c 2 denotes the gravitational radius throughout). 
With this potential, and using a constant distribution of 
specific angular momentum, l(r) = est, the effective poten- 
tial is: 



-GMbh 
R-r 



2r' 3 



dr. 



(3) 



Finally, we note that with this potential, the frequency of 
small radial oscillations for a perturbed circular orbit (i.e., 
the radial epicyclic frequency) is given by 



1 
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GM(r-3r g ) 
r(r — r g ) 3 



1/2 



(4) 



2.2 Numerical method 

For our simulations, we have us e d Smo oth Particle Hydro- 
dynamics (SPH, see iMonaehanl 1^992) for a review), in a 
two-dimensional version using cylindrical coo rdinates. We 
refer the reader to lLee fc Ramirez-Ruizl J2002T) for details of 
the implementation. There is no physical viscosity present 
in the code, only the usual artificial viscosity to model the 
presence of shocks. 

Initial conditions are generated by distributing N fluid 
elements, over the torus volume after specifying values for 
M, l(r) and <3?o, and relaxing them for several dynamical 
times in order to obtain a distribution close to equilibrium 
(see Figure 2). This configuration is then evolved in time 
with the desired perturbation to study its dynamical be- 
haviour. 

Our initial configurations are thus non-accreting tori 
fully contained within their Roche lobe. Specifically, the in- 
ner edge of the disc and the inner Lagrange point, L\ are 
located at rj n = 2.60r 9 and Tl 1 = 2.25r 9 respectively. Dur- 
ing the dynamical evolution described below, no accretion 
takes place, so the potential produced by the black hole is 
unaltered. 

The analysis of the data is carrried out by performing 
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Figure 2. An example of a thick torus in hydrodynamical equi- 
librium. This snapshot shows the relaxed initial conditions. Each 
point represents a gas cell. The locus of maximum density is at 
r = 4.25r 9 . 

the Fourier decomposition of the main hydrodynamic vari- 
ables as functions of time (e.g., the position of the center of 
the disc, the maximum and mean densities, and the various 
total energies). 



3 INTRODUCING A PERTURBATION 

The perturbation acting on the disc can be considered as 
arising from the central object, through its magnetic field, 
a deformation on its surface or a changing radiation field 
(in the case of neutron stars) or as the emission of gravi- 
tational waves from an accreting black hole. Instabilities in 
the accretion disc itself are another possible source of time 
variability, which can induce oscillations in the fluid. In ei- 
ther case, these would presumably be more intense at small 
radii. If the spin of the central object is involved in produc- 
ing the perturbation, its amplitude will vary and repeat at 
intervals given by the inverse of the spin period, AT = \ jv a . 

We have chosen, then, a perturbation which induces an 
acceleration in the disc given by: 

a per t = —r/a g ■ exp (~~~^r~) ' sm (2iri/ a t) . (5) 

Here a g is the acceleration due to gravity, r D is the outer edge 
of the torus and j; < 1 is a parameter that modulates the 
strength of the perturbation. The exponential term decays 
on a scale 5r ~ R, the radial extension of the disc, thus re- 
producing the desired behaviour for this perturbative force, 
which will be strong near the inner radius and weak in the 
outer regions. This acceleration induces radial oscillations in 
the disc, which can be Fourie r-analyzed to extract the main 
freque ncies as was done by iLee. Abramowicz fc Kluzniakl 
i2004h recently for a slender torus. 



4 RESULTS AND DISCUSSION 

A Fourier transform of the total kinetic energy in the torus 
in a typical calculation is shown in Figure 3. The black hole 
mass in this case was 2.5 solar masses, and the perturbation 
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Figure 3. FFT spectra for the kinetic (top) and total potential 
energy (bottom), for a torus orbiting a 2.5Mq black hole. From 
left to right, the first peak corresponds to the perturbation in- 
troduced at 200 Hz, the second to the radial epicyclic frequency, 
shifted to lower frequecies due to the size of the torus (the first 
acoustic p— mode), and the third to 1.5 times the second. 



frequency was fixed at v a = 200 Hz, which clearly shows up 
prominently in the spectrum as a narrow peak. 

Two additional broad features are clearly seen, centered 
at v\ ~ 300 Hz and V2 ~ 450 Hz, and are, to the limit of our 
resolution, in a 3:2 ratio. The radial epicyclic frequency for 
a test particle at the locus of maximum density (4.25r g ) is 
in this case k = 426 Hz (see equation[IJ. Trial runs with tori 
of different radial extent (Rubio -Herrera, in preparation) 
and a com parison with the work of|ZanottLRe zzolla fc Fond 
<2003ft and lRezzolla. Yoshida fc Zanottil J2003T) show clearly 
that the lower of the two peaks is simply the epicyclic fre- 
quency at the locus of maximum density, shifted to lower 
values because of the finite extent of the torus, i.e. the fun- 
damental acoustic p-mode. The second peak at higher fre- 
quency is the second in the sequence 2:3:4:... of acoustic 
p-modes. 

We draw from these results the following conclusions. 

• Global modes of o scillation in thick tori, as studi ed in 
the relativistic regime flZanotti. Rezzolla fc FonHl2003h . can 
be efficiently excited by a localized perturbation, affecting 
only the inner edge of the torus strongly. The characteristic 
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frequencies apparent in the Fourier decomposition of the 
total internal energy exhibit the sequence 2:3:... . 

• The lower of the two frequencies in the series is the 
first acoustic p-mode, closely tied to the radial epiclyclic 
frequency for test particles, and shifted to lower frequencies 
because of the finite extent of the torus. 

The nature of the excitation mechanism, as stated 
above, is left as an open question. It could be some dis- 
turbance associated with the pulsar spin frequency in the 
case of neutron star systems, or oscillations in the disc 
itself which excite these modes. In either case, the mod- 
ulation in the X-ray lightcurve will presumably occur in 
the innermost regions of the accretion flow. Since its in- 
ner boundary is defined most likely by the effects of strong 
gravity, one would expect that the frequencies would scale 
inversely with the mass of the central objec t, as is in- 
deed the case (iMcClintock fe Remillardl I2004T) for X-ray 
binaries — and may in principle allow for a mass deter- 
mination in the intermediate mass black hole candid ates 
J Abramowicz. Kluzniak. McClintock fe Remillardll2004) . 

Two potentially limiting simplifications in this study de- 
serve justification. First, in the context of LMXBs, the ratio 
of disc mass to black hole mass is very low, Md/MsH <C 1. 
Hence it is reasonable to suppose that the gravitational po- 
tential of the disc is negligible when compared to that of 
the black hole, as we have assumed. Second, we have not 
considered the effects of magnetic fields. This is simply be- 
cause we wish to carry out a study of global, purely hy- 
drodynamical modes, and does not impl y that MHD ef- 
fects are negligible or irrelevant. Recently. iKatol i2004) has 
obtained interesting results in the context of kHz QPOs 
through three-dimensional MHD calculations performed in 
a pseudo-Newtonian potential. 

Finally, one may question the choice of a constant dis- 
tribution of specific angular momentum within the torus, 
assumed here for simplicity and as a first step. How- 
ever, we note that in viscous, hydrodynamical flows, even 
if the injection of matter at large distances occurs at 
nearly Keplerian values of the angular momentum, the 
flow near the compact object may follow a flatter distribu- 
tion lllgumenshchev fe ^AjarMnowicz]^)9g|) . The Pa paloizou- 
Pringle instability ilPar^alorzou^rpTmglel Il984(l appears 
when non-axisymmetric perturbations act on non-accreting 
discs with constant angular momentum. When accretion is 
taken int o account, th is instability is suppressed, as was 
shown by [Blaej il987T> . This leads to flat distributions of 
angular momentum near the central object and a power law 
distribution in the outer region of the torus. We thus believe 
it is reasonable to assume a constant angular momentum for 
the discs as a first approximation. 

How the oscillations of the fluid in the disc may trans- 
late into variations in the X-ray lightcurve and be observed 
as kHz QPOs is ultimately still unsolved, and requires more 
complex physical processes than those included here. As- 
suming that oscillations at such frequencies do in fact occur 
due to the presence of inhomogeneities in the accretion flow, 
the luminosity modulation with account of spacetime curva- 
ture in the vicinity of a black hole has been investigated b y 
ISchnittman fc Bertschingerl l)2004h and ISchnittmanl \2004\ . 
Their results indicate that a range of variation would in fact 
be reflected in the X-rays, which encourages the investiga- 



tion of simple modes of fluid oscillation. In future work we 
will discuss the behaviour of discs with non-constant angu- 
lar momentum subject to various perturbations. 
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